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1. Recent results from HERA 

We begin by discussing new results from HERA. In autumn 2003 the ac- 
t — ' celerator started the second phase of its ep collider programme (HERA-II) . 

■ The e + p and e~p data collected by the HI and ZEUS experiments since 

then, were taken with longitudinally polarised positron and electron beams 
for the first time. Recent results from HERA-II, related mostly to the ex- 
ploration of this new feature of the collider, are discussed below. The other 
HERA result presented was an update of the high- a; measurement of the NC 
cross sections by the ZEUS collaboration 1 . In a discussion session devoted 
C~| , to the longitudinal structure function, Ft, , it was highlighted that as a test 

of different theoretical models a measurement of Ft,(x,Q 2 ) over as wide 
a range of x and Q 2 as possible would be very useful 2 , and both HI and 
ZEUS expressed a strong interest in running HERA at a low proton energy 
C$ ' for about three months, which would allow a direct Ft, measurement 3 . The 

experimental and theoretical status of the charm and beauty contributions 
to Fi arc summarised in the Heavy Flavour working group summary. 

Measurements of charged current (CC) deep inelastic scattering (DIS) 
with polarised leptons on protons allow tests of the V — A structure of CC 
interactions to be extended into the high-Q 2 regime. The polarisation de- 
pendence of the CC cross sections is fixed within the Standard Model (SM). 
Specifically, the SM predicts, from the absence of right handed charged cur- 
rents, that the CC e + p (e~p) cross section is proportional to the fraction 
of right handed positrons (left handed electrons) in the beam: 



* e cc P (V e ) = (l±V e )* e c C p (V e = 0), 



where the longitudinal polarisation V e is equal to (Nr — Nl)/(Nr + Nl) 
with Nr (Nl) being the number of right (left) handed leptons in the beam. 
HI and ZEUS published the first measurements of the CC e + p total cross 
sections for positive and negative longitudinal polarisation, and presented 5 



1 



February 2, 2008 6:33 Proceedings Trim Size: 9in x 6in 



dis06procsum 



2 



Charged Current e p Scattering 




e'p -> vX 
• H1 Data 2005 (prel.) 
: H1 Data 98-99 
' ZEUS Data 04-05(prel.) 
' ZEUS Data 98-99 



e + p -> vX 
• H1 Data 
' ZEUS Data 



20- 



40- 




-0.5 



0.5 



P. 

Figure 1. CC e+p and e~p cross sections versus V e compared to the SM prediction. 

the CC e~p total cross sections. The value of the polarisation was typically 
between 30% and 40%. The measured integrated CC cross sections, cr^ c p , 
quoted in the range Q 2 > 400 GeV 2 and y < 0.9, are shown in Fig. 1. 
The measurement of the unpolarised total cross section, in the same phase 
space, based on data collected until 2000 (HERA-I) is also shown. The 
measurements are compared to SM predictions based on the HI PDF 2000 
parametrisation 6 . The measurements agree with SM predictions and ex- 
hibit the expected linear dependence as a function of the polarisation. Lin- 
ear fits provide a good description of the data and their extrapolation to 
the point V e = —1 (V e = 1) yield a fully left (right) handed CC cross sec- 
tion for e + p (e~p) interactions which is consistent with the vanishing SM 
prediction. The corresponding upper limits on the total CC cross sections 
exclude the existence of charged currents involving right handed fermions 
mediated by a boson of mass below 180 — 208 GeV at 95% confidence level, 
assuming SM couplings and a massless right handed v e . 

The NC single differential cross sections da/dQ 2 have been measured 
at HERA-II for e + p and e~p scattering with positively and negatively po- 
larised lepton beams 7 . The SM predicts a difference in the cross section for 
leptons with different helicity states arising from the chiral structure of the 
neutral electroweak exchange. In Fig. 2 the ratio of cross sections for pos- 
itively and negatively polarised lepton beams is shown separately for e + p 
(left plots) and e~p (right plots) scattering data. In both cases the ratio is 
found to be consistent, within experimental uncertainties, with unity at low 
Q 2 , indicating little dependence of the cross section on beam polarisation. 
The normalisation uncertainties of the measurements are not included in 
the errors bars, but are indicated by the dashed lines on the upper figures. 
At higher Q 2 , the data have a tendency to deviate from unity. For e + p 
scattering the data indicate that right-handed positrons yield a larger cross 
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Figure 2. The ratios of the NC cross sections for leptons with positive and negative 
polarisations at HERA-II measured by HI (top) and ZEUS (bottom) as function of Q 2 , 
for e+p (left) and e~p (right). Solid lines correspond to the SM expectation. 



section than left-handed ones, whereas for e~p scattering the data show the 
opposite behaviour. This is consistent with the SM expectation shown as 
the solid curve in Fig. 2. 

At HERA, the NC interactions at high Q 2 receive contributions from 
7Z interference and Z° exchange. Thus, the NC data can be used to ex- 
tract the weak couplings of the light u and d quarks to the Z° boson. The 
CC cross section data help disentangle the up and down quark distributions 
and allow a determination of a propagator mass M prop from the Q 2 depen- 
dence of the cross sections. M prop measured at HERA in the space-like 
region can be compared with direct W boson mass measurements obtained 
in the time-like region by the Tevatron and LEP experiments. Combined 
electroweak and QCD fits at HERA have been performed 8 ' 9 , taking into 
account the correlation between the electroweak parameters and the par- 
ton distribution functions (PDFs). The ZEUS collaboration extend their 
NLO QCD fit 10 to inclusive NC and CC data, inclusive jet data in DIS 
and dijets in photoproduction, to determine simultaneously the PDFs and 
electroweak parameters. This fit also includes, for the first time, the e ± p 
polarised NC and CC double-differential HERA-II cross sections. The re- 
sulting PDFs are very similar to those obtained from the fit to HERA-I 
data, with slightly improved uncertainties on the u-quark at high x, due 
to additional data. In this analysis, the fits have been performed by fixing 
either the u or the d quark couplings to their SM values (fits iv-ad-PDF 
and v u -a u -PDF). The resulting one-sigma contours are shown in Fig. 3. 
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Figure 3. Contours of the 68% confidence level (CL) on the weak neutral current cou- 
plings of u (left plot) and d (right plot) quarks to the Z° boson. 

The results are consistent with the electroweak Standard Model and the 
precision is better for the u quark as expected. Comparison to the HI fits 
to unpolarised HERA-I data shows that, while the uncertainty on the axial- 
vector couplings stays about the same, the precision of the determination 
of the vector couplings is improved by a factor of 2-3 due to additional 
sensitivity of the polarised NC data. Figure 3 also shows the results of the 
HI fits in which the vector and axial-vector couplings of u and d quarks are 
fitted simultaneously (fits v u -a u -Vd-ad-PDF) and similar results obtained 
recently by the CDF experiment and at LEP. The HERA determinations 
have comparable precision to that from the Tevatron and resolve any sign 
ambiguity and the ambiguities between v u and a u of the determinations 
based on observables measured at the Z° resonance. Exploiting the Q 2 
dependence of the charged current data, the propagator mass has been 
measured to be M prop = 82. 8 ± 1.5 ± 1.3 GeV, which is in agreement with 
the direct measurements of the W boson mass. 

The amount of data collected at HERA-II is already greater than that 
of HERA-I. In particular, a significant increase of integrated luminosity 
is achieved in the e~p mode, from w 15 pb _1 to more than 100 pb _1 . 
At HERA-I the statistics of the e~p data was a limiting factor for the 
precision of the xF% determination. Profiting from the enlarged statistics 
and reduced systematic uncertainties, the previous measurement of xF$ 
has been updated 7 using HERA-II 2003-2005 e ± p NC cross section data 
at high Q 2 . Fig. 4 (left top) shows the comparison of the unpolarised 
e~p and e + p HERA-II HI data for three different Q 2 values. At high 
Q 2 , the NC cross section in e~p scattering is significantly larger than that 
in e + p scattering due to the different sign of the xF$ contribution to the 
cross section for different leptons polarities. xF$ determined by ZEUS 
from the difference of e~p and e + p cross sections is shown in Fig. 4 (right). 
The dominant contribution to xF^ arises from the "fZ interference, which 
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Figure 4. The NC reduced cross sections a NC (x,Q 2 ) measured by HI at HERA-II 
(left top), the structure functions xF% from ZEUS (right) and the averaged HI structure 
function xF^ for Q 2 = 1500 GeV 2 (left bottom) compared to QCD expectations. 



allows xF^ z to be extracted according to xF^ ~ — xF^Q 2 +M^) / (a e K,Q 2 ) 



7l z 



by neglecting the pure Z exchange contribution, which is suppressed by 
the small vector coupling v e . Here, k^ 1 = 4Af^/M|(l - M^/Mf), M z , 
Myy are the vector boson masses and a e is the axial-vector coupling of the 
electron. This structure function is non-singlet and has little dependence on 
Q 2 . The measured xF% at different Q 2 values can thus be averaged taking 
into account the small Q 2 dependence. The averaged xF^f , determined by 
HI for Q 2 = 1500 GeV 2 , is shown in Fig. 4 (left bottom) compared to the 
QCD fit result. The structure function xF$ determines both the shape 
and magnitude of the valence quark distributions independently of the sea 
quark distributions. 

2. FNAL, JLAB, RHIC and new methods for data analysis 

In this section, we summarise a number of new results from experiments at 
Fermilab, JLab and RHIC, and discuss several new approaches for deter- 
mining PDFs and their uncertainties. 

Recently, the importance of reducing the proton PDF uncertainties, es- 
pecially at high x, has become increasingly apparent, both for precision 
SM studies and new physics searches at hadron colliders. The CDF and DO 
collaborations presented a number of QCD and electroweak measurements 
that could provide additional constraints on the PDFs 11,12 . Of particular 
importance is the, poorly known, high- a; gluon distribution, which can be 
constrained using hadron collider jet data. CDF and DO presented new 
measurements of inclusive jets 11,12 in pp collisions. The data agree well 
with NLO QCD predictions over the measured Ft range, which extends to 
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Figure 5. The distribution of xs~ at Q 2 = 16 GeV 2 . 

~ 600 — 700 GcV. The precision of the data, which is limited by the jet 
energy scale uncertainty, is at a level where these measurements should pro- 
vide significant additional constraints on the high- a; gluon. In particular, in 
the most forward region measured by CDF (1.6 < \y\ < 2.1) experimental 
uncertainties are smaller than those from the PDFs. Note that measure- 
ments over a wide range of rapidity are important to discriminate between 
PDF effects and new physics signals. CDF and DO also presented a num- 
ber of other results 11 ' 12 , including measurements of the forward-backward 
charged lepton asymmetry from — ► l^v decays, which constrains the 
d/u ratio, a new measurement of the forward W — ► ev cross section and 
studies of Z — > fi/j, Z — > r e T} ll Z + jets and dijet azimuthal decorrelations. 

NLO QCD calculations of jet cross sections are CPU intensive, render- 
ing their use in QCD fits challenging. Global fit groups 13 ' 14 use LO cross 
sections + fc-factors to approximate the NLO result. At this workshop, 
an approach was presented 15 , in which an NLO QCD program is used to 
calculate a grid of weights in (x,Q 2 ), which can be convoluted with any 
PDF to give fast NLO predictions' 1 . The accuracy of the grid-computed 
cross sections is improved, without significant cost to CPU, by using co- 
ordinate transformations and high order interpolation between grid points. 
This technique will allow the rigorous inclusion of HERA, Tcvatron and 
LHC jet data in future QCD fits. 

NuTcV presented the final measurement of the difference between the 
strange and antistrangc quark distributions 17 . This is a topic of particu- 
lar interest, since it has been speculated 18 that a non-zero difference - or 
strange asymmetry - could explain the almost 3er difference between the 
NuTcV sin 2 9w result 19 and the world average. NuTeV have performed 
the first complete NLO QCD analysis of CC vN and i>N scattering, with 
two oppositely charged muons in the final state, giving direct access to the 



a A similar grid technique was also presented 16 in the Hadronic Final States session. 
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strange quark content of the nucleon. Figure 5 shows the distribution of 
xs~(x) = xs(x) — xs(x) extracted from the fit, which is positive at moder- 
ate x, such that S~ = f£ xs~(x)dx = +0.00196 ± 0.00143. NuTeV arc now 
updating their sin 2 9w result, which originally assumed S~ = 0, in light of 
this result, and also taking into account other recent measurements, such 
as the BNL K^ 3 measurement 20 . Note also, that the data prefer an asym- 
metry which forces s~(x) negative at low values of xq < 0.004, in contrast 
to some theoretical models, which suggest larger values of xq. 

The E03-103 experiment at JLab have performed a new measurement 
of the EMC effect in light nuclei 21 . The characteristic x-depcndcncc of the 
ratio <jaIo<i (i-c. the so-called "EMC effect"), has been well studied. While 
broad features can be explained, no single model has successfully described 
the effect over all x. To date, much of the experimental effort has been 
concentrated on heavy nuclei (A > 4). However, calculations of the EMC 
effect in 3 He and 4 Hc predict large differences in both magnitude and x- 
dcpcndcncc. E03-103 have measured inclusive electron scattering, using a 
5.8 GeV beam, providing the first measurement of the EMC effect in 3 Hc for 
x > 0.4, and improving upon existing measurements 22 of 4 He. The results 
indicate that the effect in 4 Hc is similar to that in carbon and the effect 
in 3 He could be larger than expected, although the strength of the latter 
conclusion is limited by large uncertainties in proton excess corrections. 

JLab also presented results on the phenomenon of parton-hadron duality 
in semi-inclusive scattering. Duality, for example the fact that the average 
behaviour of nuclear resonances exhibits the scaling behaviour expected in 
pQCD, is well established in inclusive scattering, but has yet to be studied 
to the same extent in semi-inclusive data. The merit of semi-inclusive pro- 
cesses lies in the ability to identify individual quark species in the nucleon 
by tagging specific mesons in the final state. The E00-108 experiment at 
JLab presented 23 preliminary measurements of semi-inclusive pion electro- 
production, eN — > e'n ± N, using a 5.5 GeV electron beam on proton and 
deuteron targets, for 1.8 < Q 2 < 6.0 GeV 2 , 0.3 < x < 0.55, and elasticities, 
z = E-x/v, in the range 0.35 — 1. The results indicate that the ratio of the 
unfavoured to favoured fragmentation functions, ^qr = ^ N +jJj Nl ' 1 > shows 
no dependence on x and a smooth slope as a function of z, in accordance 
with the expectations of the onset of duality. It was argued that these data 
are also consistent with the idea that pion electroproduction is related to 
low energy factorisation. If so, future semi-inclusive measurements could 
give excellent access to the flavour structure of the nucleon. 



February 2, 2008 6:33 Proceedings Trim Size: 9in x 6in 



dis06procsum 



Vs m =200 GeV 



p+p->7T°+h ± +X d+Au^7T°+h ± +X 



• n° (<7)>=4.00) 
□ h" (»)=3.2) 
o h" (?)=2.2) 

fan 



'_a 



Ti mesons <7?> = 4.00 

shadowing (KKP) 

- — shadowing(Kretzer) 
■ • • multiple scattering 



.5 0.2 

TO 
D 



Q_ 

0)0.2 




S=0.100±0.014 



25<E.<30 GeV 



2.5 5 
S=0.154±0.024 



0.2 - 



-S=0.020±0.013 



25<E„<30 GeV 



3 2.5 5 
: S=0.093±0.040 



30<E„<55 GeV 





p T (GeV/c) <P« ~ <Plcp 

Figure 6. Left: nuclear modification factor, i?dAu = 2 x 197 ^d+Au/op+p i as a function 
of pt for mesons at (rj) = 4.00 and charged hadrons at smaller pseudorapiditics. 
Error bars are statistical, while the shaded boxes are point-to-point systematics. Right: 
azimuthal correlations between forward 7r° and charged hadrons at mid-rapidity. 

The STAR collaboration have searched for parton saturation effects at 
RHIC 24 . The term saturation describes the idea that at very small x, 
the gluon density is so large, that gluons begin to recombine, leading to 
a slowing of the growth and giving rise to non-linear terms in the evolu- 
tion equations. The state of saturation is often referred to as the Colour 
Glass Condensate (CGC). Nuclear environments, such as RHIC, are ideal 
places to search for signs of saturation, due to the A 1 / 3 enhancement of 
the parton densities in nuclei, relative to the proton. STAR presented re- 
sults on forward 7r° production in p+p and d+Au collisions, in the range 
3-0 < i] < 4.2, at = 200 GeV. Measurements at forward rapidities 
ensure that low values of gluon- a; are probed. The p+p yields generally 
agree with NLO pQCD predictions, while the d+Au yield is suppressed at 
forward rapidities, as shown in Fig. 6 (left). It was argued that the pr 
dependence of the d+Au yield is consistent with a model treating the Au 
nucleus as a CGC for forward pion production 24 . The azimuthal corre- 
lations of the forward ir° with charged hadrons at mid-rapidities (Fig. 6 
(right)) show a recoil peak in p+p that is suppressed in d+Au collisions 
at low E n ; in qualitative agreement with a gluon saturation picture of the 
Au nucleus. Note, however, that the PHENIX experiment at RHIC do not 
see significant differences in the p+p and d+Au azimuthal correlations 25 . 
STAR will soon make further measurements, utilising a new forward meson 
spectrometer, that could elucidate the source of the observed suppression. 

Several presentations were made in which alternative approaches to PDF 
determination were explored. The NNPDF collaboration presented the 
latest results on a neural network (NN) approach to PDF fitting 26 . This 
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approach attempts to avoid any bias introduced from a choice of functional 
form for the PDFs and should provide more faithful estimates of the PDF 
uncertainties. The NNPDF best-fit of the non-singlet structure function, 
F 2 JVS (x, Q 2 ) = F%(x, Q 2 ) - F 2 d (x, Q 2 ), agrees with the predictions of other 
PDF sets, and gives a better description at large x, as well as larger error 
bands in the region where there is no data. Work is now in progress to 
construct a full set of PDFs. A related approach was presented by the 
SOMPDF group 27 . This technique is based on a specific class of NN, called 
the Self Organising Map (SOM). It was argued that SOMs allow better 
control of systematic bias, by replacing the fully automated procedure of 
standard NNs, with an interactive fitting procedure. Finally, the prospects 
for determining PDFs and their uncertainties using a Bayesian statistical 
approach, was discussed 28 . This work is at an early stage, but it will be 
extremely interesting to see the results when they become available. 

3. Theory 

The main topics were updates in determinations of proton and nuclear par- 
tons distributions and the ever increasing sophistication (complication?) 
of theory. The former included the implementation of new heavy flavour 
prescriptions and/or NNLO corrections and new data in fits. The latter 
included recent developments in small- a; resummations, with hopefully the 
beginning of detailed understanding and phenomenology. Matters are get- 
ting more involved in the very small x region with non-linear effects, and 
the saturation scale seems to be becoming ever smaller. 

Starting with fits to partons distributions, Guffanti presented a NNLO 
analysis of non-singlet parton distributions 29 u v (x,Q 2 ) and d v (x,Q 2 ) by 
fitting to F%' d (,x,Q 2 ), for x > 0.3 and Ff - F 2 d . The result is a s (M|) = 
0.1134 ± 0.0020, in good agreement with other NNLO fits to only structure 
function data. Alekhin made his fits more global 30 by including E605 Drell- 
Yan data 31 and E866 Drell-Yan ratio data 32 . The fit has no problems with 
this data, and it improves the accuracy on the high- a; sea and gives his 
first real constraint on u — d. In fact using A^ 2 = 1 the sea uncertainties 
are very small. The theoretical input to the DIS fit is the massless NNLO 
QCD corrections for the light quarks and gluons (splitting 33 and coefficient 
functions 34 ' 35 ; account of the heavy quarks up to C>(a|) 36 in FFNS; account 
of the target-mass corrections, Fermi-motion in deuterium, and twist 4 
terms; and the massless 0(a|) corrections 37 to the coefficient functions. 
There is disagreement with the definition of NNLO regarding heavy quarks 
- this is conventionally only NLO corrections in the FFNS scheme. The last 
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step is only part of a full N 3 LO correction - lacking the splitting functions 
and not necessarily indicative except at high x. Reasonable stability is 
claimed down to Q 2 = 0.5GeV 2 for 0.06 < x < 0.12 - perhaps possible 
since most corrections are small in this region. 

There are some issues with new data in the fits. Most interesting was 
the NuTeV structure function data 38 which is larger than the similar CCFR 
data 39 at high x and is useful for flavour separation. Comparison to this 
data relies on nuclear corrections. A determination of these was reported by 
Kumano 40 . CTEQ find these NuTeV data difficult to fit 41 , whilst MRST 
only obtain a good fit by severely reducing the nuclear correction 42 , im- 
plying that this is different for neutrinos than for charged leptons. The 
important information in the data is for the region x < 0.3 which is not too 
sensitive to the nuclear corrections, but the problem caused much interest. 
Recent CHORUS data 43 are more similar to CCFR data. There are also 
some changes in global fit procedures. CTEQ 41 include all HERA data 
and fit directly to cross-sections for first time which requires Fl(x,Q 2 ) 
at high y. They also implement a new heavy flavour prescription 44 (see 
the Heavy Flavour session summary). It was noted that the new fit over- 
shoots the high-?/ data, with the match between theory and data requiring 
the use of the photoproduction background systematic error at more than 
2 cr, i.e. use of the systematic errors can remove the high-?/ Fl(x,Q 2 ) 
turnover. MRST implement a full NNLO VFNS 45 and use the Drell-Yan 
cross-sections 35 for the first time, leading to (provisional) full NNLO par- 
tons with uncertainties 42 . The improvement in the charm procedure af- 
fects the gluon compared to the approximate MRST2004 NNLO fit 14 . The 
quality of the fit at NNLO is consistently better than NLO, but not for 
Drell-Yan data. There is a definite tendency for as(M|) to go down with 
all the changes - at NLO a s (M 2 ) = 0.121 and at NNLO a s {M 2 z ) = 0.119. 

There were presentations on various ways to include small- a; resum- 
mations from the BFKL equation on top of the fixed order expansion. 
White showed a resummation of ln(l/a;) terms along with running cou- 
pling corrections 46 , which produces mainly analytic results with small nu- 
merical corrections. The procedure includes the quark-gluon splitting P qg 
and the full implementation of a heavy flavour VFNS 47 and gives full LO 
rcsummcd partons. Colferai outlined an approach which is also applicable 
for processes with two hard scales. It too includes the running coupling but 
also resummation of collincar singularities: 

Jdk 2 (k 2 )-i- 1 K n (k 2 ) = X n h), X n (7) ~ l/7 2 " + \ l/(l-7) 2 " +1 
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The evolution variable is s/(QQo) with conjugate variable N. Consid- 
eration of changes of evolution variable to s/Q 2 (DIS) and s/Q 2 , lead 
to the rcsummation 48 X n(i) ~ 1/(7 + N/2) n+1 , 1/(1 - 7 + N/2) n+1 . 
The natural calculations arc in DIS scheme with the incoming gluon off- 
shell - k 2 = Qq 7^ 0. Colferai looked at the transformation to the MS 
scheme 49 , showing that the effect is not large. Forte presented an approach 
based on: duality - one has x{l(N, as), as) = N, "f(x(os, M), as) = N 
i.e. the Q 2 evolution and x evolution are dual 50 (perhaps not the most 
important issue); explicit imposition of momentum conservation; inclu- 
sion of the running of the coupling 51 ; and symmctrisation, i.e. letting 
1/M n 1/(M + N/2) n , 1/(1 - M) n — > 1/(1 — M + N/2) n . Overall this 
leads to a resummcd P gg (x, Q 2 ) 52 . These methods are all rather different in 
details, but now have many similarities, and produce similar results for P gg 
at NLO in the resummation, with a dip below the LO splitting function for 
x ~ 0.001, and a slow growth not setting in until x ~ 10~ 5 . White also ex- 
amines phenomenology, at LO in the rcsummation only - the impact factor 
required for P qg is not yet known at NLO. The analysis gives a better fit 47 
than NLO-in-as in terms of dF^{x,Q 2 )/d\a.Q' 2 , but the enhancement of 
the evolution is too great at small x and the resulting gluon and Fl(x, Q 2 ) 
are too small at moderate x, suggesting that we need the full NLO gener- 
alization. Colferai 49 and White 53 both examine the improvement to P qg . 
The two approaches are qualitatively similar but contain different higher 
order information: White estimates NLO corrections to the impact factor 54 ; 
Colferai has a rcsummation of P gg beyond NLO via the collinear resumma- 
tion. Both suggest the effects of the NLO resummation beyond fixed order 
in the quark sector are small but significant, as seen in Fig. 7. 




Figure 7. The splitting function P qg in two approaches including small x rcsummations. 

There was an update on non-linear corrections at very small x, with 
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various discussions on how to include the saturation corrections, and more. 
Soyez considered the nonlinear evolution equation in rapidity Y extended 
to include fluctuations as well as recombination 55 , i.e. the equation is 

9r(TW) = olMbfkl ® (T (n) ) - aM BFKL <g> (T {n+1 ^) + aa 2 s K g> (T^" 1 )), 

where T m is the amplitude for m dipole scattering. Hatta demonstrated 
the origin of the fluctuations via a formal derivation of a Bremsstrahlung 
Hamiltonian 56 which gives the evolution of n-dipole densities. Soyez showed 
that the inclusion of the fluctuations leads to a dispersion about the satu- 
ration scale 57 . One can be sensitive to saturation effects when nominally 
far away from the saturation scale, or conversely, when naively near this 
scale, fluctuations to even lower x have no effect since saturation is already 
reached, while those to higher scales reduce saturation effects. The net 
effect is to move the onset of saturation in the dipole cross-section consid- 
erably downwards, as seen in Fig. 8, which compares the dipole cross-section 
with and without fluctuations. Kutak modifies the LO BFKL equation to 
include a non-linear recombination term 58 and also investigates impact pa- 
rameter dependence 59 . Moreover, high- a; effects in the gluon evolution are 
accounted for (though not in the gluon-quark impact factor, i.e. P qg )- The 
impact parameter dependence reduces the effect of the nonlinear term and 
lowers the saturation scale in Q 2 for fixed x, i.e. peripheral collisions see 
smaller densities. The conclusion from some rough phenomenology is that 
we "hardly see the effect of saturation" at HERA. 




-8 -6 -4 -2 2 4 

log(l/r 2 ) 

Figure 8. Comparison of the dipole cross-section from a fit with only saturation cor- 
rections (solid) and one which also allows for fluctuations (dots). 

To summarise, there is little agreement in global fit analyses. Not every- 
one wants to go to NNLO, and not everyone agrees how to do it in detail. 
However, it seems we are now at the stage where NNLO parton analyses 
are essentially complete and reliable and should be done. They work a little 
better than NLO in general. There are rather similar results coming from 
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all groups working on small- a; resummations to be used in conjunction with 
fixed order calculations, though differences in procedures. The effect of the 
resummations is moderate until very small x. Empirically, their inclusion 
can improve the fit a little over NLO. There is progress in nonlinear small- x 
equations, e.g. fluctuations, which always seems to be pushing the satura- 
tion scale lower. It would be good if this approach could match on to higher 
x better - at present it is usually confined to some unknown small x region, 
and is missing higher x corrections. Overall there is lots of improvement in 
how to calculate using different techniques. However, there is not enough 
idea yet where each approach is applicable/needed. We still need better 
(real) phenomenology and, of course, more useful data. 
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